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Abstract 

We consider the Potts model and the related bond, site, and mixed site-bond percolation prob- 
lems on triangular-type and kagome-type lattices, and derive closed-form expressions for the critical 
frontier. For triangular- type lattices the critical frontier is known, usually derived from a duality 
consideration in conjunction with the assumption of a unique transition. Our analysis, however, 
is rigorous and based on an established result without the need of a uniqueness assumption, thus 
firmly establishing all derived results. For kagome-type lattices the exact critical frontier is not 
known. We derive a closed-form expression for the Potts critical frontier by making use of a ho- 
mogeneity assumption. The closed-form expression is new, and we apply it to a host of problems 
including site, bond, and mixed site-bond percolation on various lattices. It yields exact thresholds 
for site percolation on kagome, martini, and other lattices, and is highly accurate numerically in 
other applications when compared to numerical determination. 

PACS numbers: 05.50.+q, 02.50.-r, 64.60. Cn 
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I. INTRODUCTION 



An outstanding problem in lattice statistics is the determination of the critical frontier, 
or the loci of critical point, of lattice models. Of special interest is the g-state Potts model 
[l| and its associated lattice models j^. For q = 2 it is the Ising model, and for q = 1 the 
Potts model generates the percolation problem Q| including the bond 4], site {5], ane mixed 
site-bond percolation. However, except for the simple square, triangular and honeycomb 
lattices [?] and some special lattices essentially of a triangular-type [8], the determination 
of the Potts critical frontier in general has proven to be elusive. 





FIG. 1: (a) Triangular- type lattice, (b) Kagome-type lattice. Shaded triangles possess Boltzmann 
weights (P). 



In this paper we consider the Potts model on two general classes of lattices, the triangular- 
and kagome-type lattices shown in Fig. [TJ Shaded triangles in Fig. [T] denote interactions 
involving 3 Potts spins r 1; r 2 , r 3 = 1, 2, q with the Boltzmann weights 



W A {1, 2, 3) = A + B(S 12 + 5 23 + 5 31 ) + CS 123 , 
W v (l, 2, 3) = A' + B'(5 12 + 5 23 + 5 31 ) + C'5 l23 



(1) 



where 5^ = S^fa, Tj), 5i 23 = 5 12 S 23 S 3 i, and A, B, C, A', B', C are constants. In ([ID, we have 
assumed interactions isotropic in the 3 directions of a triangle. The extension of our analysis 
to anisotropic interactions is straightforward and will not be given. Special cases of shaded 
triangles are the "stack-of-triangl e" , or subnet, lattices shown in Figs. [2] and [3] that have 

9|, [lO, llll, ll2|- We refer to these stack-of-triangle lattices as subnet 



been of recent interest 

lattices. The lxl subnet lattices are the triangular and kagome lattices themselves. We 
shall call a kagome-type lattice with mxm down-pointing and nxn up-pointing subnets 
an (m x m) : (n x n) subnet lattice, or simply an [m x m) : (n x n) lattice. Examples of 
these kagome-type subnet lattices are shown in Fig. [3j 




(a) 



(b) 



(c) 



FIG. 2: Triangular subnet lattices, (a) lxl lattice (triangular), (b) 2x2 lattice, (c) 3 x 3 lattice. 




FIG. 3: Kagome subnet lattices, (a) (1 x 1) : (1 x 1) lattice (kagome). (b) (2 x 2) : (2 x 2) lattice, 
(c) (3 x 3) : (3 x 3) lattice, (d) (1 X 1) : (2 x 2) lattice. 

Partition functions for the two lattices in Fig. [1] are 

i 

Z tli (q;A,B,C) = ^l[W A (i,j,k), ( 2 ) 

n=l A 
Q 

Z^fa A, B,C; = Y,[Y[W*{i,j,k)] ■ [l[W v (i',j',k')], (3) 



Ti = l A 



V 



where the products are taken over the respective up- and down-pointing shaded triangles. 
The critical frontier of the triangular-type lattice of Fig. Ufa) has been known from earlier 



works [1 



14 



151 ]. but the critical frontier for the kagome-type lattice of Fig. QJb) is open. 



For q = 2, we can replace Potts spins r by Ising spins a = ±1 and shaded triangles by 
triangles (lxl subnet) with an Ising interaction Kj (= K/2). To determine Kj, we write 



ij 



^(l+OiOj), 



Wa(1,2,3) = exp[2K I (6 12 + 6 23 + S 31 )}. 



Equating (j4]) with Wa in ([I]), one obtains after a little algebra 



(4) 



e 4K ' = (A + 3B + C)/(A+ B). 



(5) 
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It follows that the partition functions ([2]) and ([3]) are completely equivalent to those of the 
triangular and kagome Ising model. 

For q > 3, the shaded triangles can be replaced by any triangular network having 2 
independent parameters. An example is the mapping shown in Fig. [T2] in Sec. IIII Gl 

Parameters A, B, C for a given Potts subnet can be readily worked out. For the lxl 
triangle, for example, one has 

W A (1, 2, 3) = exp[K(5 12 + 5 23 + 5 31 )} (6) 

from which one obtains 

A = 1, B = v, C = 3v 2 + v 3 , (triangle), (7) 

where v — e K — 1. For the 2x2 subnet, one obtains in a similar fashion 

A = 3v 5 + 21w 4 + (50 + Aq)v 3 + 33qv 2 + 9q 2 v + q 3 , 
B = v 7 + 7v 6 + 22v 5 + (30 + 2q)v i + l®qv 3 + gV, 

C = v 9 + 9v 8 + 33v 7 + 63v 6 + (54 + 3q)v 5 + 9qv\ (2x2 subnet). (8) 

Expressions of A, B, C for 3x3 and 4x4 subnets are derived and given in a subsequent 
paper [if]], hereafter referred to as II. 

The structure of this paper is as follows: In Sec. [Til we consider the triangular subnet 
lattices and apply the rigorously known critical frontier to various models including mixed 
site-bond percolation. In Sec. Ill we consider the kagome-type lattice, and derive a closed- 
form expression for its critical frontier on the basis of a homogeneity assumption. We show 
that this critical frontier is exact for site percolation on the kagome, martini, and other 
lattices, and is highly accurate in other applications. The accuracy of the critical frontier 
will be closely examined in paper II. 



II. TRIANGULAR-TYPE LATTICES 



In this section we consider triangular- type lattices of Fig. [U^a). 

The Potts model on the triangular-type lattice was first studied by Baxter, Temperley 
and Ashley [l3] in the context of a Potts model with 2- and 3-site interactions. Using a 
Bethe-ansatz result on a 20-vertex model on the triangular lattice due to Kelland 



13, 



they showed that the partition function (|2J) is self-dual, and derived its self-dual point which, 
in the language of the interaction (JTJ), reads 

qA = C. (9) 

This self-dual trajectory was later re-derived graphically by Wu and Lin [ijj]. However, as 
is common in duality arguments, an additional assumption of a unique transition is needed 
to ascertain that (19]) is indeed the actual critical frontier. 



However, Wu and Zia 



established subsequently in a rigorous analysis that (jSJ) is 



indeed the critical frontier in the 'ferromagnetic' regime 

2B + C>0, 3B + C>0. (10) 

It can be verified that the condition (1101) holds for (j7j) and (|S]), so the critical frontier qA = C 
is exact. Applications of (jUJ) to the martini and other lattices have been reported in [8j. The 
duality relation of the triangular Potts model with 2- and 3-spin interactions has also be 
studied by Chayes and Lei [l9(] with several rigorous theorems on the phase transition proven. 



A. Ising model 

In Sec. I, we have established that for q = 2 any triangular-type lattice is reducible to 
a triangular Ising lattice with interaction Kj given by (jSJ). Indeed, using (jSJ), the known 
critical point e AKl = 3 of the triangular Ising model reduces to the critical frontier 2A = C 
as expected. 

For the Ising model on lxl and 2x2 subnet lattices with interaction Kj, we set 
q = 2, v = e 2Kl — 1 in (J7|) and (JHJ), and obtain from 2 A = C the critical point 

x c = y/3, lxl subnet (triangular lattice) 

= V5, 2x2 subnet, (11) 

where x = e 2Kl . Using expressions of A and C given in II for 3x3 and 4x4 subnets, 
we obtain similarly 

x 8 - 5x 6 - x 4 - 19x 2 -8 = 0, x c = ^5 + ^+^(50 + 18^) 

= 2.404 689 372, 3 x 3 subnet 
x 12 - 5x w - x 8 - 22x 6 - 53a; 4 - 125a; 2 - 51 = 0, x c = 2.467 648 033, 4 x 4 subnet. (12) 
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B. Bond percolation 



It is well-known that bond percolation is realized in the q = 1 limit of the Potts model 
under the mapping v = pj (1— p), where p is the bond occupation probability 0, 4]. Therefore 
the percolation threshold is given simply by C = A. Thus using (J7J) for A and C for the 



triangular lattice, one obtains the well-known [20 



21 



221 ] critical frontier for bond percolation 



1 - 3p + p 3 = 0, or p c = 2 sin(7r/18) = 0.347 296 355. (triangular lattice). (13) 

For the 2x2 subnet lattice we use © and obtain 

1 - 3p 2 - 9p 3 + 3p 4 + 45p 5 - 72p 6 + 45p 7 - 12p 8 + p 9 = 0, 

or p c = 0.471 628 788 (2 x 2 subnet lattice). (14) 

In a similar fashion using expressions of A and C given in II, we obtain 

1 - 3p 3 - 18p 4 - 39p 5 + 77/ + 309p 7 - 198p 8 - 1406p 9 + 315p 10 + 9303p n - 23083p 12 
+28707p 13 - 22047p 14 + 10959p 15 - 3462p 16 + 636p 17 - 52p 18 = 0, 

p c = 0.509 077 792 (3 x 3 subnet lattice), (15) 
1 - 3p A - 30p 5 - 114/ - Q3p 7 + 636p 8 + 1940p 9 + 741p 10 - 14283p n - 26541p 12 
+78759J9 13 + 189279J9 14 - 370589p 15 - 1229877p 16 + 2829339p 17 + 6938691j9 18 - 41655363j9 19 
+96750306p 20 - 143421123p 21 + 152405700p 22 - 121438416p 23 + 73822093p 24 - 34270647p 25 
+ 11994555p 26 - 3073478p 27 + 545409p 28 - 60012p 29 + 3089p 30 = 0, 

p c = 0.524 364 822 (4 x 4 subnet lattice). (16) 



These findings agree with those of Haji-Akbari and Ziff [12J deduced from a duality con- 
sideration. As aforementioned, our derivation now ascertains that these thresholds are the 
exact transition points. 



C. Potts model 

The exact critical threshold for the Potts model on triangular- type lattices is 1^, or 
qA = C . Using expressions of A and C given in ([7]), one obtains the known critical frontier 

3v 2 + v 3 = q, (Potts model on triangular lattice). (17) 



For the 2x2 subnet lattice one uses (jSj) and obtains the critical frontier 

v 9 + 9v 8 + 33v 7 + 63v 6 + 54t» 5 - 12qv 4 - (50g + Aq 2 )v 3 - 33gV - 9q 3 v - q 4 = 0. (18) 

Solutions of (TTTT) and (|T8|) and those of the 3x3 and 4x4 subnet lattices are tabulated in 
Table[T]for q — 1, 2, 3, 4, 10. Note that the q = 1 solutions are related to the bond percolation 
thresholds (IT3]) -(ITfi|) by e A ' c = 1/(1 -p c ). 



TABLE I: Exact Potts threshold e c for triangular-type subnet lattices. 





q = l 


q = 2 (Ising) 


q = 3 


q = A 


q = 10 


Triangular lattice 


1.532 088 885 


V3 


1.879 385 241 


2 


2.492 033 301 


2x2 


1.892 608 790 


V5 


2.493 123 120 


2.706 275 430 


3.602 637 947 


3x3 


2.036 982 609 


2.404 689 372 


2.674 398 828 


2.895 419 068 


3.808 005 450 


4x4 


2.102 451 724 


2.467 648 033 


2.731 876 784 


2.946 645 097 


3.820 754 228 



D. Site percolation 

n 

Kunz and the present author [5j have shown that site percolation can be formulated as a 
q = 1 limit of a Potts model with multi-site interactions. The Kunz-Wu scheme is to consider 
a reference lattice with multi-spin interactions, and regard faces of multi-spin interactions 
as sites of a new lattice on which the site percolation is defined. The critical frontier of the 
Potts model on the reference lattice then produces the site percolation threshold for the new 
lattice. This scheme of formulation can be extended to mixed site-bond percolation. 

1. Site percolation on the triangular lattice: 

Consider as a reference lattice the triangular lattice with pure 3-site interactions M 
marked by dots shown in the left panel of Fig. |H The dots form a triangular lattice shown 
in the right. The Kunz-Wu scheme now solves the site percolation on the triangular lattice. 
We have 

W A (l,2,3) = e MS ™ = l + m5 123 , (19) 

or A = 1, B = 0, C — m = e M — 1. Writing m = s/(l — s), where s = 1 — e~ M is the 
site occupation probability, the exact critical frontier A = C now yields immediately the 



well-known site percolation threshold [20 



21 



22 



23j for the triangular lattice, 



1/2. (20) 



FIG. 4: Site percolation on the triangular lattice. 



2. Site percolation on kagome lattice: 

Consider the reference 2x2 subnet lattice with pure 3-site interactions M denoted by 
dots shown in the left panel of Fig. [5j The Kunz-Wu scheme then maps the reference Potts 
model into site percolation on the kagome lattice as indicated in the right. 

Now for a 2 x 2 subnet containing 3 dots as in Fig. [5l we have 



A = q 3 + 3qm, B = m 2 , C = m 3 



(21) 



where m = e M — 1. Writing m = s/(l — s) and setting q — 1, the rigorous critical frontier 
A = C yields the critical condition 1 — 3s 2 + s 3 = 0, leading to the known exact result 



20. 



2l| 



s kag = x _ 2s in(7r/18) = 0.652 703 644. 



(22) 





FIG. 5: Site percolation on the kagome lattice. 
3. Site percolation on (1 x 1) : (n x n) lattices: 

The above scheme of mapping can be extended to site percolation on (1 x 1) : (n x ri) 
lattices for general n. The example of Fig. [5]is n = 1, and the n = 2 lattice is shown in Fig. 
|6j The reference lattice (not shown) for n = 2 consists of 3 x 3 subnets with 

A = q 7 + Qq 5 m + 15g 3 m 2 + (lAq + 3q 2 )m 3 + 3m 4 , 
B = q 2 m 3 + 2(q + l)m 4 + m 5 , 

C = 3m 5 + m 6 , (23) 
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where m = e M — 1, M is the 3-site interaction. After setting q = 1 and m = s/(l — s), the 
critical frontier qA = C becomes 

(l-2s 2 )(l + 2s 2 -3s 3 + s 4 ) = 0, (24) 

yielding the exact threshold 

s c = l/y/2, (1 x 1) : (2 x 2) kagome site percolation. (25) 
The exact critical threshold for higher (lxl) : (n x n) lattices can be similarly worked out. 



FIG. 6: Site percolation on the (1 x 1) : (2 x 2) lattice. 



III. KAGOME-TYPE LATTICES 

We consider in this section the case of the kagome- type lattices of Fig. Hfb). 

The critical frontier of the Potts model on kagome-type lattices has proven to be highly 
elusive. On the basis of a homogeneity assumption, however, the present author [7|] has 
advanced a conjecture on the critic al p oint for the kagome lattice. The conjecture has 



since been closely examined [24j, |25j, |26|, |27| and found to be extremely accurate. Here we 
extend the homogeneity assumption to general kagome-type Potts lattices. For continuity 
of reading, we first state our result in Sec. IIII A\ and present the derivation in Sec. IIII Gl 

A. A closed-form critical frontier and homogeneity assumption 

For Potts model on kagome-type lattices described by the partition function (j3J), the 
critical frontier under a homogeneity assumption is 

{q 2 A + 3qB + C)(q 2 A' + 3qB' + C) - 3{qB + C)(qB' + C) - (q - 2)CC' = 0. (26) 

Remarks: 



1. Despite its appearance, the critical frontier fl26l) actually contains only 3 independent 
parameters (Cf. fl55|) below). 

2. The expression (l2"6l) is exact for g = 2. 



B. Ising model 

We first show that ( l26i) is exact for q = 2. We have already established that the partition 
function ([3]) is precisely that of the kagome Ising model. For completeness we now verify 
that the critical frontier (j2"6"|) also gives the known kagome critical point. 

For symmetric weights A = A', B = B', C = C, the kagome Ising model has a uniform 
interaction Kj and the critical point is known to be at e iKl = 3 + 2\/3~ {28]. It is readily 
verified that, by using ([5]) for Kj, this critical point gives rise to precisely the q = 2 critical 
frontier (1261) . namely, 

AA + 6B + C = V3(2B + C). (27) 

The proof can be extended to the kagome-type model with asymmetric weights. 

Critical thresholds of kagome-type Ising subnet lattices computed from (127]) are tabulated 
in Table IIHI For the kagome and (2 x 2) : (2 x 2) Ising lattices, for example, we use ([7]) 
and ([8]) for A, B, C with q = 1, v = x — 1, x — e 2Kl . This gives 



x 4 - 6x 2 - 3 = 0, x c = -\/3 + 2V3 

= 2.542 459 756, (kagome lattice) , 



x s - 8x 6 - 6x 4 - 32s 2 - 83 = 0, x c = y 2 + V3 + \J 12 + lOv 7 ^ 

= 3.024 382 957, (2 x 2) : (2 x 2) lattice. (28) 



C. Bond percolation 

For bond percolation threshold on kagome-type subnet lattices, we again use (|2T)|) with 
the substitution of q = 1 and v = p/(l — p), where p is the bond occupation probability. 
Using (J7|) and (jHJ), we obtain 

1 — 3p 2 — 6p 3 + 12p 4 - 6p 5 +p 6 = 0, p c = 0.524 429 717 (kagome lattice), (29) 
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TABLE II: Bond percolation threshold p c for (m x m) : (n x n) lattices for m, n < 4. 



1 J CX) U L 1 V V. 


TVi i k wnvl^ 

A i.J.10 VV Ul XV 


TViimprip^l HptPTTniria tion 

J. > U.J.J.i.V'J. J.V_>OiJ. uvl J. J. J. J. 1. J. CI Ij J. Ull 


Kagome 
(1 x 1) : (2 x 2) 
(1 x 1) • (3 x ,3 s ) 
(1 x 1) : (4 x 4) 


0.524 429 717 
0.570 882 620 
599 798 340 
0.592 017 120 


0.524 404 99(2) [27] 
0.570 866 51(33) [29] 


(2 x 2 s ! ■ (2 x 21 
(2 x 2) : (3 x 3) 
(2 x 2) : (4 x 4) 


600 870 248 
0.610 916 740 
0.614 703 624 


600 862 4(1 0) M 1 1 


(3 x 3) : (3 x 3) 
(3 x 3) : (4 x 4) 


0.619 333 485 
0.622 473 191 


0.619 329 6(10) [11] 


(4 x 4) : (4 x 4) 


0.625 364 661 


0.625 365 (3) [11] 



1 - 3p 3 - 12p 4 - 12p 5 + 63p 6 + 60p 7 - 330/ + 423p 9 - 264p 10 + 8Ap n - lip 12 = 0, 

p c = 0.570 882 620, (1 x 1) : (2 x 2) lattice, (30) 

1 - 3p A - 18p 5 - 39p 6 + 30p 7 + 273/ + 264p 9 - 1785p 10 - 126p n + 8232p 12 - 16236p 13 

+ 16359p 14 - 9948p 15 + 3708p 16 - 786p 17 + 73p 18 = 0, 
p c = 0.600 870 248, (2 x 2) : (2 x 2) lattice. (31) 

Bond percolation thresholds computed from (126]) for (m x m) : (nxn) lattices are tabulated 
in Table ILT1 We also include in Table ILT1 numerical determinations of p c for the (1 x 1) : (2 x 2) 



29J and (n x n) x (n x n), n = 2, 3, 4 lj| lattices by Ziff and Gu using simulations, and 



of the kagome lattice by Feng, Deng and Blote [27] from a transfer matrix analysis. The 
comparison shows that (1261) is accurate to within one part in 10 5 . 



D. Potts model 



Critical thresholds for the Potts model on kagome-type subnet lattices computed from 
( |26l) are tabulated in Table [TTT1 For the kagome lattice itself, for example, we have A = A' = 
1, B — B' — v, C — C — 3v 2 + v 3 , v = e K — 1, and ( 1261) gives the critical frontier 

v 6 + Qv 5 + 9v 4 - 2qv 3 - 12qv 2 - 6q 2 v - q 3 = 0, (kagome lattice). (32) 
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TABLE III: Potts threshold e Kc for ka gome-type subnet lattices. 



1 .aft i ( 


(7—1 

q — i 


ff / I G1TKT 1 

q — z ^ismgj 


q — o 


n — 4 
q — i 


n — 1 n 


kagome lattice 
(1 x 1) : (2 x 2) 
(~\ y 1V fty 1) 

l±/\ J- 1 . 1 O /\ O J 

(1 x 1) : (4 x 4) 


2.102 738 619 
2.330 364 713 
9 4Q8 74D 9fiD 

2.451 083 242 


2.542 459 757 
2.821 281 889 
9 903 973 669 
2.928 442 860 


2.876 269 226 
3.186 678 923 
3 960 483 758 
3.276 998 285 


3.155 842 236 
3.489 096 458 
3 553 390 863 
3.562 314 883 


4.355 385 241 
4.761 529 399 
4 764 908 41 
4.739 553 252 


f9 x 2\ ■ (9 x 91 
(2 x 2) ■ f3 x 3) 
(2 x 2) : (4 x 4) 


9 505 450 QOQ 
2 570 143 984 
2.595 404 635 


3 094 389 957 
3 082 166 484 
3.098 624 716 


3 481 055 307 
3 454 087 416 
3.378 293 046 


3 71 7 691 699 
3 757 519 846 

1 U 1 "7 -L L7~ V7 

3.761 399 505 


5 01 6 339 590 
5 004 155 712 
4.984 524 206 


(3 x 3) : (3 x 3) 
(3 x 3) : (4 x 4) 


2.626 971 274 
2.648 818 511 


3.133 002 727 
3.147 204 863 


3.497 087 416 
3.416 364 328 


3.712 498 867 
3.796 037 357 


4.992 841 134 
4.973 931 010 


(4 x 4) : (4 x 4) 


2.669 262 336 


3.160 721 132 


3.598 289 910 


3.639 241 821 


4.954 642 401 



The critical frontier (|32|) for the kagome lattice was first obtained by the present author some 
30 years ago [?], 3(3] by using the homogeneity assumption described in Sec. IIII Gl Compar- 



ison of the thresholds computed from (I3"2"|) for q = 1, 3, 4 with Monte Carlo renormalization 
group findings has shown that the accuracy of (1261) is within one part in 10 5 241 ]. 



E. Site and site-bond percolation 




FIG. 7: Site percolation on the 3-12 lattice. 
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The 3-12 lattice is the lattice shown in the right panel of Fig. [7J To formulate site 
percolation on the 3-12 lattice, we consider the reference (2 x 2) : (2 x 2) lattice with pure 
3-site interactions shown in the left. Let the 3-site interactions of the up- and down-triangles 
be, respectively, Mi and M 2 . One finds 



A = q i + 3qm 1) B = m\, C = m\, 



A' = q 3 + 3qm 2 , B' 



m 



2" 



C = m 3 



2- 



(33) 



where m 1 = e Ml — 1, m 2 = e M% — 1. Setting q — 1, m% — s\/ (1 — si), m 2 = s 2 /(l — s 2 ), where 
Si and s 2 are the respective site occupation probabilities for the 3-12 lattice, the critical 
frontier f l26j) gives 



1 - 3(sis 2 ) 2 + (sis 2 ) 



0. 



(3 — 12 site percolation). 



For Si = s 2 = s, this yields the known 20|, |23j critical frontier 1 — 3s 4 + s 6 = 0, or 

y/l - 2sin(7r/18) = 0.807 900 076. 



,3-12 



(34) 



(35) 



Using the relation s^ 9 = (s^ _12 ) 2 23], we have therefore derived the exact kagome and 3-12 
site percolation thresholds, and demonstrated that ( l26i) is exact in this instance. 
2. Site percolation on the martini lattice. 





FIG. 8: Site percolation on the martini lattice. 



The martini lattice [3jJ is the lattice shown in the right panel of Fig. [HJ To generate a site 
percolation on the martini lattice, we start from the (lxl): (2x2) reference lattice with 
3-site interactions shown in the left. Denote the 3-site interactions of up- and down-pointing 



triangular faces by Mi and M 2 respectively and write m\ = e Ml — 1, m 2 = e m — 1. We have 



,M 2 



A 
A* 



q 3 + 3qmi, B = ml, C = m\ 
1, B' = 0, C = m 2 . 



(36) 
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Setting q = 1, mi = si/(l — si), m 2 = S2/(l — s 2 ), with si and s 2 the respective site 
occupation probabilities, fl26|) gives the critical frontier 



l-{3sl-sl)s 2 
This is a known exact result 



0, (martini site percolation). 



(37) 



26 



32 



33j . and is another example that the critical frontier 
(|26|) is exact. For s 2 = 1, the percolation reduces to that on the kagome lattice, and (!37|) 
gives the threshold (1221 . For uniform occupation probability S\ = s 2 = s, fl37|) becomes 
1 - 3s 3 + s 4 = and gives the exact solution s ™ artini = 0.764 826 486. 
3. Site-bond percolation on the honeycomb lattice. 

No exact result is known for the site and site-bond percolation on the honeycomb lattice. 
Owing to the intrinsic interest of a percolation process on a simple Bravais lattice, the 
problem of honeycomb site percolation has attracted considerable attention for many years. 
There now exists a host of highly precise numerical estimates on the threshold for site 



percolation on the honeycomb lattice 



11 



23 



27] . 









p 






^ 2 > 




FIG. 9: Site-bond percolation on the honeycomb lattice. 

Consider the more general mixed site-bond percolation on the honeycomb lattice with 
site occupation probabilities S\ and s 2 and bond occupation probability p shown in the 
right panel of Fig. [9j The relevant reference lattice can be taken as shown in the left with 
edge-interactions K and 3-site interactions Mi and M 2 . To make use of ( |26|) . we adopt the 
scheme of devising up- and down-pointing triangles as indicated in Fig. [TlTb) below. This 
gives 



A = (q + v ) 3 + (q + 3v )m 2 , B = w 2 m 2 , C = v 3 m 2 , 
A' = 1, B'= 0, C = mi, 



(38) 



where v — e K — 1, rrij = e * — 1, i — 1, 2. Setting q — 1, v — p/(l — p), rrii — Sj/ (1 — Sj), we 
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obtain from (126]) the critical frontier for the mixed site-bond percolation as 
(3p 2 — p 3 )siS2 = 1, (honeycomb site — bond percolation). 



(39) 



When Si = S2 = 1, (139]) is exact since it gives the known honeycomb bond percolation 



threshold 1 - 3p 2 + p 3 = [20J 
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22l | . When p — 1, f[39|) gives the threshold 

sis 2 = 1/2 



(40) 



which is exact for s 2 = 1, as the site percolation reduces to one on the triangle lattice with 
the critical point (l2"0l s„ = 1/2. But for st = a, = a, (llOl) gives s c = l/y/2 = 0.707 106 781 
differing from accurate numerical estimates of s c = 0.697 040 2 [27{| and s c = 0.697 041 3 
13J . The critical frontier (126]) is therefore a close approximation in this instance. 

The site-bond percolation has also been studied by simulations by Ziff and Gu for s\ = 
S2 [ll] and for p = 1 Q|. Their results indicate fl39|) works better for site occupation 
probabilities ~ 1. 

4- Site-bond percolation on the kagome lattice. 





FIG. 10: Site-bond percolation on the kagome lattice. 

Consider the mixed site-bond percolation on the kagome lattice with site and bond oc- 
cupation probabilities s and p shown in the right panel of Fig. [101 The reference lattice is 
shown in the left having edge interaction K and 3-site interaction M. Regard the reference 
lattice as a kagome-type with the partition function ([3]). One has 



A 
B 
C 
.4' 



q{q 2 + 3m) (g + v f + 3m 2 (g + v f 
m 2 (q + v ) 2 + m 3 v, 
m 3 {3v 2 + v 3 ), 

1, B' = v, C' = 3v 2 + v 3 , 



m 



(41) 
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where v = e K — l,m = e M — 1. Substituting (14TI) into (126]) and setting q = l,v = p/(l — 
p), m = s/(l — s), one obtains the critical frontier 

1 + 3s 2 (l - 3p + 2p 3 - p 4 ) + s 3 (-3 + 9p - 3p 2 - 12p 3 + 15p 4 - 6p 5 + p 6 ) = 0, 

(kagome site — bond percolation). (42) 

For p — 1, (1421) becomes 1 — 3s 2 + s 3 = which gives the exact critical threshold (1221) for 
the kagome site percolation. For s = 1, (j4"2l becomes 1 — 3p 2 — 6p 3 + 12p 4 — 6p 5 + p 6 = 0, 
or p c = 0.524 429 717, in agreement with (1291) . 

F. The 3-12 lattice 




(a) (b) 



FIG. 11: (a) The 3-12 Potts lattice, (b) The 3-12 lattice as an asymmetric kagome-type lattice. 

The 3-12 lattice is the lattice shown in Fig. [TTT a) with interactions K, K\, K 2 . To make 
use of (126]) we regard the lattice as one of the kagome-type consisting of large up-pointing 
triangles (dotted lines) and small down-pointing triangles as indicated in Fig. fTTY b) . Then 
we have (see also Eq. (5) of |8|]) 

A = (q + v) 3 + 3{q + v)(q + 2v)v 2 + 3(q + 3v)vj + (q + 3v)v 3 2 

B = v 2 [{q + v)v 2 + 3v 2 2 + v 3 } 

C = v 3 {3v 2 2 +v 3 ) 

A' = 1, B' — v h C = 3v\ + vf, (43) 

where v = e K — l,vi = e Kl — l,v 2 = e K2 — 1. Substituting (|43j) into (|26|) . we obtain the 
critical frontier (re-arranged in a symmetric form) 

(q+v) 3 (h 1 +3qv 1 +q 2 )(h 2 +3qv 2 +q 2 )-3(qv 2 +v 3 )(h 1 +qv 1 )(h 2 +qv 2 )-(q-2)v 3 h 1 h 2 = 0, (44) 

where hi = 3vf + vf, i — 1, 2. 
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For the 3-12 Ising model with uniform interactions Kj, we set q = 2, V\ = v-i = v 
e 2Kl — 1, and (pElj) simplifies to 



(V3-l)v 2 -2u-4 = 0, 



(45) 



yielding the known exact critical point e 2K ' = |(3 + v 7 ^) + \ (6 + 5^/2 = 5.073 446 135 



in agreement with Utiyama 



34| and Syozi [28 1 



For Potts model on the 3-12 lattice with uniform interaction K, ()44l) gives 

v 9 + 6v 8 + 3(3 - q)v 7 - g(32 + q)v 6 - q(75 + 30q)v 5 

- g 2 (lll + 12g)w 4 - 2g 3 (41 + q)v 3 - 36q 4 v 2 - 9q 5 v - q e = 



(46) 



This gives the critical point 



v + 1 



3.852 426 158, q = 1 

5.073 446 135, q = 2 (exact Ising result) 

6.033 022 515, q = 3 

6.857 394 828, q = 4. 



(47) 



The accuracy of the prediction (1471) will be examined in paper II. 

For bond percolation on the 3-12 lattice we set q = 1 and write v = p/(l — p),V\ = 
Pi/(1 — f 2 —Vil (P2—P2), where p,Pi,P2 are the respective bond occupation probabilities. 
Then (j4"4"j) gives the critical frontier 

1 - 3p 2 ( Pl +pj -pD( P 2 +pl -Pi) + P 3 (3p 2 - 2pl)(3p 2 2 - 2p\) = 0, 

(3 — 12 bond percolation). (48) 

This expression has been conjectured recently by Scullard and Ziff as a non-rigorous 
extension of the exact bond percolation threshold of the martini lattice. In the uniform case 
Pi = P2 = P, (SSI) becomes 



1 -p + p +p - Ip +4f> = 0, p c 



0.740 423 317, 



(49) 



which is also given by p c = 1 — e Kc using e Kc for q = 1 in ( j47j ) . Compared to the numerical 
determination of_p c = 0.740 421 95(80) by Ziff and Gu [ll| and the value p c = 0.740 420 81 
by Parviainen 35j, the accuracy of the homogeneity determination (149!) is seen to be well 
within one part in 10 5 . 
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For the mixed site-bond percolation on the 3-12 lattice, it is temptin g to use the kagome 



23j. This gives the 



critical frontier (H2l and replace s by s 2 p as argued by Suding and Ziff 
critical frontier 

1 + 3s 4 (j» - 3p 2 + 2p 4 - p 5 ) + s 6 (-3j9 3 + 9p 4 - 3p 5 - I2p 6 + I5p 7 - 6p 8 + p 9 ) = 0. (50) 

For p = 1 the pure site percolation, this becomes 1 — 3s 2 + s 6 = which is exact. For s = 1 
the pure bond percolation, (1501) gives 

1 + 3p 2 - 12p 3 + 9p 4 + 3j9 5 - 15/ + 15p 7 - 6p s + p 9 = (51) 

with the solution p c = 0.747 882 617. The small difference between fj4*9l and fl5T|) reflects 
the approximate nature of the kagome site-bond critical frontier 



G. Critical frontier and homogeneity assumption 

We now derive the critical conjecture (1261) using a homogeneity assumption. 

In the partition function (T5]), we replace the two Boltzmann weights by 

q 

W&(1, 2, 3) = Fx e L(5 ii' +5 22'+' 5 33') e ^<V 2 /3' j 

<1 

W v (l,2,3) = F' x Y e i,(5 "' +<522 ' +5 33') e ^i'2'3' ; (52) 

s l' s 2' s 3 = 1 

as indicated graphically in Fig. [12j Equating ( 1521) with (TjQ), we find 

A = Fx [(q + if + (q + 3i)n], 
B = F xfn 

C = Fx fn, (53) 

where I = e L — 1, n = e N — 1 and similar relations for A', B', C with {F, £, n} — > {F', £', n'}, 

£' = e L ' -l,n' = e N ' -1. 

Solving (|53|) for £, n, F, one obtains 

B' 

(qB + Cf 



n 



AC 2 — qB 3 — 3B 2 C 
B 3 (AC 2 - qB 2 - W 2 C) 



F C 2 (qB* + Cf ' 
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FIG. 12: Graphical representation of (f52|) . 



and similarly one obtains i',n',F' in terms of A',B',C. The kagome-type lattice now 
becomes the one shown in Fig. [131(a). 





(b) 



FIG. 13: (a) The kagome-type lattice after transformation (|52j) . (b) The lattice dual to (a) 



The duality relation of Potts models with multi-site interactions has been formulated by 
Essam 22] (see also [2]). Following Essam, the dual to the lattice in Fig. ITBla) is the one 
shown in Fig. [TBT b) with 



e K = (l + q/£)(l + q/£'), 
e M = l + q 2 /n, 



1 + (f/ri 



(55) 



where the interaction K is the dual to the two interactions L and V in series. We therefore 
are led to consider the Potts model on the triangular lattice shown in Fig. [131( b) . where M 
and M' are 3-site interactions. 

For M' = 0, the partition function is Z tri (q; A, B, C) given by ([2]) with 



W A (l,2,3) = e 



K (Sl2+523 +6 3 l) p M8l23 
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or 

A=l, B = e K -l, C = e 3K+M -3e K + 2. (56) 
The exact critical frontier in this case is known. It is qA = C, or 

e 3K+M _ 3e K + 2 = g ^ (5?) 

For M' 7^ the critical frontier is not known. However, the critical frontier must be symmet- 
ric in M and M' . We now make a homogeneity assumption requiring M and M' to appear 
homogeneously in the exponent of (1571) . The simplest way to do this is to extend (1571) to 

e 3K+M+M' _ 3e K + 2 = q (5g) 

The substitution of expressions of K, M and M' in ( 155]) and (1531) into ( 1571) now leads to (l2"6l) . 
IV. SUMMARY 

We have considered the g-state Potts model and the related bond, site, and mixed site- 
bond percolation for triangular- and kagome-type lattices. For triangular-type lattices we 
obtained its exact critical frontier in the form of (Q without the usual assumption of a 
unique transition. We then applied the exact critical frontier in various applications. For 
kagome-type lattices we obtained a new critical frontier (1261) by making use of a homogene- 
ity assumption. We established that the new critical frontier is exact for q = 2 and for site 
percolation on the kagome, martini, and other lattices. For the Potts and bond percolation 
models for which there is no exact solution, the new critical frontier gives numerical values 
of critical thresholds accurate to the order of 10 -5 . For mixed site-bond percolation, the ho- 
mogeneity assumption gives rise to critical frontiers which are accurate when site occupation 
probabilities are ~ 1. 

In summary, we emphasize that applications of the critical frontiers (El) and ( 1261) are not 
limited to those reported in this paper. They can be extended to numerous other lattice 
models having a triangular or kagome symmetry, and thus they open the door to a host of 
previously unsolved problems. 
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